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We present a theory for recent experiments on the quantum spin liquid candidate material 
Na4lr308- In particular, we provide an explanation of large enhancement and high sensitivity of the 
Wilson ratio and magnetic susceptibility across the metal-insulator transition (MIT). We consider 
an extended Hubbard model on a hyperkagome lattice, which incorporates atomic spin-orbit cou- 
pling (SOC) and multi-orbital interactions of Ir 5d electrons. We analyze the model near MIT and 
study thermodynamic properties across the MIT. We propose that the system on the insulating side 
is a U(l) quantum spin liquid with spinon Fermi surfaces. We further show that, the magnetic sus- 
ceptibility experiences multiple enhancements from the strong SOC, reduced bandwidth and onsite 
spin-orbital exchange, while the heat capacity does not. Understanding of the currently available 
experiments and suggestions for future experiments are discussed. 
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Quantum spin liquid (QSL) is an exotic state of matter 
with fractionalized excitations and emergent gauge fields. 
QSLs on the theoretical side are well established, but it 
has been very challenging to identify the direct evidence 
for the existence of QSLs in real materials. Recently, 
however, several materials have been proposed as promis- 
ing candidates of QSL lj. While no magnetic ordering 
has been observed in these systems, there exist ther- 
modynamic signatures indicating the presence of gapless 
spin-carrying excitations with a finite charge excitation 
gap. One natural explanation of these spin-carrying ex- 
citations would be spinons in QSL phases. The central 
question in the context of QSLs, therefore, is what kind 
of spin liquid phases may be realized in these systems. 

In this work, we present a theory of a three-dimensional 
QSL that may be realized in Na^^Os, an Ir-based hy- 
perkagome lattice system. We focus on recent experi- 
ments of anomalously enhanced Wilson ratio across the 
metal-insulator transition (MIT) in this material. Using 
a semi-realistic model of Ir ti g orbitals, it is shown that 
the anomalously large Wilson ratio is a natural conse- 
quence of the QSL phase arising from strong spin-orbit 
coupling (SOC) and multi-orbital interactions. Com- 
bined with other spin liquid phenomenology, such a large 
Wilson ratio would be a fingerprint of a spin-orbital QSL, 
where the spinons have both spin and orbital characters. 

In the early experiment on polycrystalline NaziIraOs 
samples in Ref. ® the magnetic susceptibility and heat 
capacity did not show any signature of magnetic order- 
ing down to the lowest measurable temperature. NMR 
Knight shift measurement further confirms the absence 
of magnetic ordering down to 2K 3j. With the Curie- 
Weiss tempeature Qcw — — 650K, the frustration pa- 
rameter (defined as / = e ^ v with Tjy the ordering 
temperature) is greater than 300. Moreover, in the zero 
temperature limit, the magnetic susceptibility saturates 
to a finite and large constant while the heat capac- 



ity has a linear temperature dependence with a rather 
small coefficient 7 (7 = ^f). This leads to an anoma- 

lously large Wilson ratio W = \ x/ /'f at low tem- 
perature, that equals to 35 ! Most of other promis- 
ing QSL candidate materials have a Wilson ratio of or- 
der of unity at low temperatures pQ. In particular, the 
organic materials such as k-(BEDT-TTF) 2 Cu 2 (CN) 3 [I] 
and EtMe3Sb[Pd(dmit) 2 ] 2 [5] also show constant mag- 
netic susceptibility and linear-in-temperature heat capac- 
ity, often taken as an evidence for spinon Fermi surface. 
These materials are close to a MIT and it has been sug- 
gested that a QSL phase of spinon Fermi surfaces may 
arise upon gapping out charge excitations via a second 
order transition starting from the metallic side O [7]. In 
this case, the spinon Fermi surface would be a remnant of 
the electron Fermi surface and the Wilson ratio would not 
change much across MIT. While Na^^Os may also be 
close to a MIT due to the extended nature Ir c£-orbitals, a 
simple extension of this argument would not explain the 
large Wilson ratio. Clearly, a different or an additional 
physics is at work for Na 4 Ir 3 8 . 

Recently there have been new experiments on single 
crystal samples [3 . Depending on the preparation, the 
single crystal samples can be both insulating and metal- 
lic. Comparing the most conducting metallic sample and 
the insulating polycrystalline sample, one finds that, the 
heat capacity only has a very small change with 7 var- 
ied by 17% while the magnetic susceptibility is increased 
by a factor of 7. In general, for the single crystal sam- 
ples, the magnetic susceptibility is much more sensitive 
to MIT than the heat capacity, and is strongly enhanced 
on the insulating side. 

There have been several theoretical proposals about 
the nature of the disordered state on the insulating 
side (SHU]. It should be pointed out that all of the exist- 
ing proposals have assumed a spin- 1/2 Hcisenberg model 
on a hyperkagome lattice. Ref. [TUI suggested a U(l) 
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QSL with spinon Fermi surfaces based on a projective 
wavefunction study. Such a state is certainly qualita- 
tively consistent with the Fermi-liquid-like phenomenol- 
ogy in Na4lr 3 0g. The Wilson ratio in this proposal, how- 
ever, should still be order of unity because of the spin- 
rotational invariance. Ref. ED introduced a Z2 QSL with 
spinon pairing. They attribute the large Wislon ratio 
to the suppression of the heat capacity by the spinon 
pairing. Even though the singlet spinon pairing also sup- 
presses the magnetic susceptibility, they argue that this 
can be avoided if SOC energy scale is greater than the 
pairing strength. This state, if relevant for Na4lr 3 0s, 
may indicate emergent superconductivity for the metal- 
lic samples at similar temperature scales. If there is no 
superconductivity, the metallic samples are expected to 
have a much larger heat capacity than the insulating 
ones. Neither is observed in the recent experiments [3J. 

Although the ground state of the Heisenberg model on 
the hyperkagome lattice itself is an interesting theoreti- 
cal problem, the applicability of this model to the actual 
material can only be justified in the strong Mott regime 
with certain assumptions [12 . In fact, both the polycrys- 
talline sample and single crystal samples are proximate 
to a MIT. We also notice that all other QSL candidate 
materials can be modelled by either a single-band Hub- 
bard model in an intermediate or the strong Mott regime. 
On the other hand, in Na-JraOs , all the three ti g orbitals 
are involved and hence the electron interactions are or- 
bitally dependent. Moreover, since Ir is a heavy element, 
a large SOC is naturally expected|12j. An important 
question is which aspects of these ingredients give rise to 
qualitatively different thermodynamic behaviors. 

In this paper, we analyze the extended Hubbard model 
on the Ir-based hyperkagome lattice, including all the 
three ta g orbitals, SOC, and the electron interactions on 
these orbitals. We assume that the MIT is controlled 
by the correlation. We separate the multi-orbital inter- 
actions into the Hubbard-C/ interaction for the charge 
sector and the onsite spin-orbital exchange for the spin- 
orbital sector. We study the MIT in this model by a 
slave-rotor mean-field theory. The main result of this 
paper is that, the strong SOC and the correlation both 
suppresses the electron bandwidth, which effectively en- 
hances the onsite spin-orbital exchange in the insulating 
side. SOC breaks the spin-rotational symmetry and en- 
hances the bare electron/spinon magnetic susceptibility, 
and the magnetic susceptibility is further enhanced by 
the "enhanced" onsite spin-orbital exchange. While the 
spin suceptibility gets multiple enhancements, the heat 
capacity is only sensitive to the density of state on the 
Fermi surface which does not experience strong enhance- 
ment. This leads to a large enhancement in the Wilson 
ratio across the MIT. 

To write down our model, we make use of a few 
well-known facts about the microscopies related to 
Na4lr30s|iT2] The local IrOg crystal field first splits the 



5d electron states into a low-lying ta s triplet and a much 
higher e g doublet which will play no role. The atomic 
SOC acts on the three t2 g orbitals. Moreover, there may 
also exist further crystal field splitting among the three 
ta g orbitals that arises from the small distortion of the 
IrOg octahedron 2J . Including the hopping of electrons, 
the non-interacting part of the Hamiltonian is given by 



(ij) ,mn,OL 

Here d ima describes the electron with orbital m — 
xy, yz, xz and spin a J, at site i, and L (er) is the or- 
bital angular momentum (spin Pauli matrix). In Eq. ([I]), 
the single-ion anisotropy Di m captures the distortion of 
the IrOg octahedra. For the hopping, we include both 
direct and oxygen p-orbital mediated indirect electron 
tunneling. The details about different terms in Ho are 
further discussed in the supplementary material|13j. 

For the interaction, we adopt the standard multi- 
orbital interaction that includes intra-orbital repulsion 
([/), inter-orbital repulsion (U'), Hund's coupling (J) and 
pairing hopping (J'), 



U' 



(2) 



Here, n ? 



ima d ima and n i 
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the atomic limit, these four Kanamori's parameters in 
Eq. ^ have the relation U = U' + J + J' , J = J' , which 
is assumed in the following discussion. Such a multi- 
orbital interaction is certainly very complicated to deal 
with. To manifest the physics, we follow the treatment in 
Ref. 1141 and decompose this multi-orbital interaction into 
the charge part H c and the spin-orbital part H ex with 
Hi = H c + Hex and 

U 



j 



(3) 



z,m^m' ,(7(7' 



Here, we assume the average electron occupation per site 
is 5, appropriate for the Ir 4+ ion with 5d 5 electron con- 
figuration. We neglect an unimportant constant in H c . 
H c is the usual Hubbard-t/ interaction and describes 
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the energy cost for the electron charge fluctuation. T~L ex 
describes how the electrons arranges themselves among 
different spin-orbital states, i.e. onsite spin-orbital ex- 
change interaction. Since H c and H ex describe two dif- 
ferent physics, we will treat them separately. 

The Hubbard U is the largest energy scale in the inter- 
action and is treated in the first place. We first consider 
the Hamiltonian % = Hq + Ti. c and analyze the phase 
diagram within the slave-rotor mean- field treatment |15j. 
In the slave-rotor formalism, we decompose the electron 
operator into a bosonic charge rotor e l9i and a fermionic 
spinon fi ma (that carries spin and orbital quantum num- 
bers), i.e. d ima = e~ %6i Uma- With this decomposition, 
the physical Hilbert space is enlarged and we need to im- 
pose a constraint L t = J2 ma flmJima ~ 5 to S et back to 
the physical Hilbert space. Here, L.; is the angular mo- 
mentum operator conjugate with 9i, [O^Lj] = iSij. The 
Hamiltonian H expressed in the rotor and spinon vari- 
ables is further decomposed into two mean-field Hamil- 
tonians for the rotors and the spinons, respectively, 



Wr=|D £ » 3 + fcL *)+«-E e<(fl '" 9i) & 

1 (ij) 

tij,mnfi ma fj na h ^ ] fi ma fim<j- (5) 

{ij),mn,a i,m,a 

Here, h is the Langrange multiplier introduced to 
implement the constraint on average, and the mean 
field parameters Q r and Qf are given by Q r = 

E m „ CT %,mn</L CT /j„ CT )/ and Q f = (eW*-W) r , where 
the subindices indicate the mean field ground state to 
evaluate over. Here we have made a uniform mean field 
approximation such that Qf and Q r are uniform on all 
the bonds and the ground state is smoothly connected to 
the non-interacting case. When the rotor is condensed 



with 



^ 0, the spinon binds with the charge ro- 



tor, form an electron with a finite quasi-particle weight 
Z^O, and the system is in a Fermi liquid phase. When 
the rotor is uncondensed with (e lSi ) = 0, the spinons are 
deconfined, form spinon Fermi surfaces, and the system is 
in a U(l) QSL phase. We solve this mean-field Hamilto- 
nians self-consistely. The rotor condensation is signaled 
when the lowest rotor mode turns gapless. 

In the phase diagram depicted in Fig.[T] there are three 
different energy scales: SOC, correlation and bandwidth. 
SOC narrows the electron bandwidth, which effectively 
enhances the correlation effect and drives a MIT at a 
reduced critical interaction. This is the key underlying 
reason for the presence of strong correlation physics in 
5d electron systems jTHj which were long believed to be 
weakly correlated. Moreover, the correlation also sup- 
presses the bandwidth and thus enhances SOC effect. 

Using Ioffe and Larkin's relation[T5], we have for the 




FIG. 1. (Color online) The slave-rotor mean-field phase 
diagram[l7] with the Hubbard-f/ in the charge sector only. 
i CT is the a-bonding hopping amplitude Q3]. Different curves 
corresponds to different single-ion anisotropy parameter D. 
From left to right, D = 0.8^,0.4^,0.2^,0. 



heat capacity C v = C v f+C vr with C v f (C vr ) the spinon's 
(rotor's) contribution. Since the rotor's contribution is 
subdominant at low temperatures, C vr cx T 3 in the Fermi 
liquid phase, C vr oc e~ A / T (A is the charge gap) in 
the Mott insulator, we then have C v « C v f. The mag- 
netic susceptibility comes only from the spinon's contri- 
bution, i.e. X — Xf- I n the slave-rotor formalism for 
the Hamiltonian %, the spinons are essentially treated as 
free fermions. In Fig. [2j we have computed the thermo- 
dynamic quantities for two different SOCs (A = 0.8^, t a ) 
with D = and D = 0.2t a . The specific heat is 
not monotonic because the spinon/electron Fermi surface 
changes due to the presence of the anisotropy and the 
SOC. The specific heat with D = only varies slightly 
near the MIT, which is consistent with the experiments. 
Taking t a = 0.64eV and A = t a = 0.64eVQJ], we find 
7 ~ 3.63mJ/molK 2 , only slight larger than the experi- 
mental value («3mJ/molK 2 )[3]. The specific heat with 
D = 0.2t a is strongly suppressed near the MIT in con- 
trast to the experiments. This suggests the actual ma- 
terial may not have a strong anistropy due to the IrOg 
distortion. In Fig. |2j we only find a small enhancement of 
magnetic susceptibility, which is inconsistent with the ex- 
periments. Since the Wilson ratio is sensitive to magnetic 
fluctuations, and moreover, it is the magnetic suscepti- 
bility that varies significantly in the experiments, we nat- 
urally turn to the onsite spin-orbital exchange Hex that 
has not been so far included in our analysis. 

We follow Ref. HU and place the onsite exchange H ex 
into the spinon mean-field Hamiltonian % / and treat the 
Hex by a weak coupling analysis. This treatment is jus- 
tified by the fact that T-L ex describes the spin-orbital ex- 
change rather than the charge fluctuation. Moreover, to 
further assist this treatment, we also take the exchange 
coupling to be much weaker than the Hubbard U inter- 
action. The new spinon mean-field Hamiltonian Wfn is 
given by Hf n = Hf + H ex with % ex expressed in terms 
of spinon operators. Since the ground state does not 
have magnetic order, the inclusion of the exchange into 
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FIG. 2. (Color online) The specific heat, susceptibility and 
Wilson ratio obtained from slave-rotor mean-field theory with 
the Hubbard-!/ in the charge sector only. In the left (right) 
three plots, we have A = 0.8t CT (A = t a ). In all the plots, 
the red circular (blue squared) dots are data points with D = 
(D = 0.2i CT ). Dashed curves connecting the dots are for 
guidance. The arrows indicate the location for metal-insulator 
transition. The results are obtained for a finite system with 
40 x 40 x 40 unit cells at T = 0.0005i CT . 



the spinon Hamiltonian does not modify the mean-field 
phase diagram in a weak coupling approach. Moreover, 
the density of state on the Fermi surface is not modified 
by the inclusion of the exchange, so the heat capacity 
stays invariant as the case without the exchange. 
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FIG. 3. (Color online) Magnetic susceptibilities and Wilson 
ratios for A = t a , D = 0, including the spin-orbital exchange. 
In both plots, from top to bottom, J = 0.3i CT , Q.2t a , 0.1t CT , 0. 
Dashed curves connecting the dots are for guidance. The ar- 
rows indicate the location for metal-insulator transition. The 
results are obtained for a finite system with 40 x 40 x 40 unit 
cells at T — 0.0005t CT . 

The magnetic susceptibility can be strongly enhanced 
by this onsite spin-orbital exchange. Because the spin- 
rotational symmetry is explicitly broken by the SOC and 
the orbital angular momentum also contributes to the 
magnetization, the magnetic susceptibility is not just 



sensitive to density of state at the Fermi surface but is 
sensitive to the whole many-body states. This is quite 
different from a normal Fermi liquid with spin-rotation 
invariance. The onsite spin-orbital exchange favors the 
electrons to occupy different orbitals and thus increase 
the orbital angular momentum. Note that putting elec- 
trons onto a single orbital state tends to quench the or- 
bital angular momentum. The orbital angular momen- 
tum further couples to the spin via SOC. The presence 
of spin-orbital exchange strongly modifies the system's 
response to external magnetic fields. 

The resulting magnetic susceptibilities and Wilson ra- 
tios with different spin-orbital exchange are plotted in 
Fig. [3j The magnetic susceptibility is only slightly en- 
hanced by the exchange on the Fermi liquid side, while 
it is strongly enhanced on the QSL side. Here two rele- 
vant energy scales are the electron bandwidth and the 
exchange. In the Fermi liquid phase, the bandwidth 
is broad and a weak exchange does not cause much 
change. In the QSL phase, the bandwidth is suppressed 
by the correlation and SOC, which effectively enhances 
the onsite exchange. The enhanced onsite spin-orbital 
exchange, together with the SOC, leads to a large en- 
hancement of magnetic susceptibilities and Wilson ratios. 

Here we discuss additional experimental probes to de- 
tect the QSL in Na/JrsOg- It is known that Ir absorbs 
neutron which prevents neutron scattering measurement 
for small size samples. A recently developed experimen- 
tal technique, resonant inelastic X-ray scattering, can 
overcome this difficulty in iridates and may provide an 
opportunity to observe the spinon continuum in the low- 
energy spin sector [20]. Transport measurement such as 
thermal conductivity may also provide an indirect evi- 
dence for the spinon Fermi surfaces if a metal like behav- 
ior is observed pi]. 

To summarize, we propose that the ground state in the 
insulating Na4lr30g is a U(l) quantum spin liquid with 
spinon Fermi surfaces. We show that the anomalously 
large enhancement and high sensitivitiy of the Wilson 
ratio across the metal-insulator transition is due to the 
combined effects of strong SOC, correlation and onsite 
spin-orbital exchange. 
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torted and generates a single-ion anisotropy in the Hamil- 
tonian Hq. At each Ir site, there exists a 2-fold rotational 
axis and the symmetry-allowed single-ion anisotropy has 
several independent terms. We pick up a term that is 
likely to be the dominant single-ion anisotropy. For in- 
stance, the anisotropy at sublattice 4 in Fig. [4] is given 
by T,mn,a D [( lZ ) 2 ]mndl ma d in(r , where l mn = L mn and 
D > 0. Anisotropics on other sublattices can be readily 
obtained by symmetry. 

Our tight-binding model contains three parameters: er- 
bonding t a and 7r-bonding t n for direct electron tunnel- 
ing, and t it i hopping for the indirect electron tunneling 




FIG. 4. (Color online) Ir based hyperkagome lattice and its 
parent pyrochlore lattice. "1,2,- • • ,12" label the 12 sublat- 
tices. Solid bonds (in red) connect Ir sites. 



through the intermediate oxygen p-orbitals. We neglect 
the higher-order processes involving the tunneling be- 
tween two different oxygen p-orbitals. For instance, the 
hopping Hamiltonian between sublattice 4 and 5 is of the 
form, 



SUPPLEMENTARY MATERIAL 



Non-interacting Hamiltonian 

As shown in Fig. |4j the Ir hyperkagome lattice has 
12 sublattices in the unit cell. Each IrOg is slightly dis- 
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In this paper, we choose t^jt a — 0.1 and Ud/t a = 0.6 
that are close to the ones used in the band structure 
calculation in Ref.19 of the main text. 



